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ABSTRACT 



For very slow white dwarf accretors in CV's Townsley and Bildsten (2004) 
found a relation between the accretion rate M and the central temperature of 
the white dwarf Tp. According to this relation for M less than 10~^^MQyr~^ Tc 
is much lower than 10'^ K. Motivated by this study we follow the thermonuclear 
runaway on massive white dwarfs (M^^ = 1.25 — 1.40 Mq) with Tc lower 
than 10^ K, accreting matter of solar composition. We demonstrate that 
in that range of the relevant parameter space {Tc,Mwd and M) the slope 
of the relation between the peak temperatures achieved during the runaway 
and Tc becomes much steeper than its value for Tc above 10^ K. The peak 
temperatures we derive can lead to nuclear breakout from the conventional "hot 
carbon-nitrogen-oxygen" cycle. When breakout conditions are achieved the 
heavy element abundances can show a much wider variety than what is possible 
with the common enrichment mechanisms. 



Subject headings: binaries:close-novae,cataclysmic variables, nucleosynthesis, 
abundances - stars:white dwarfs 
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1. Introduction 

Classical nova explosions are a consequence of the unstable ignition of hydrogen in 
accreted envelopes on white dwarfs in close binary systems. Peak temperatures in the range 
~ 200 — 350 million K achieved in the ensuing thermonuclear runaway yield burning of 
hydrogen to helium by means of the hot CNO cycles, with the associated nuclear energy 
release ~ 5 x 10^^ erg g~^ powering the outburst. One dimensional numerical simulations 
Starrfield et.al 1974, Starrfield et.al 1985, MacDonald 1980, MacDonald 1983, Nariai et. 
al. 1980, Prialnik et.al 1978, Prialnik 1986, Iben 1982, Truran 1982, Jose & Hernanz 1998 
have established that many features of classical novae in outburst can be understood on 
this model, for representative choices of white dwarf mass {Mwd ~ ^o) and accretion rate 
(M a; 10"^° — lO~^M07/r~^ ). More details can be found in a review by Starrfield 2002. 

Observational studies of novae reveal that, while the composition of the accreted matter 
is consistent with Solar abundances, the ejecta are characterized by high concentrations 
of heavy elements (e.g. C, N, O, Ne, Mg) relative to solar abundances Gehrz et.al 1998. 
For most nova models, numerical simulations predict that peak temperatures achieved in 
the runaways do not exceed the critical value {Tcrit of ~ 4 x 10^ K) at which "breakout" 
from the hot CNO cycles can occur Wiescher et.al 1986, Yaron et.al. 2005. This has led to 
the conclusion accepted by most nova theorists that the observed enhancements reflect the 
effects of " dredge-up" of CO- or ONeMg-rich matter from the underlying white dwarf star 
Livio & Truran 1990. We show here that in a restricted, and therefore rare, region of the 
relevant parameter space ( Tc, Mwd and M ) breakout burning can lead to heavy element 
enrichment caused only by breakout burning. 

In a recent exploration of the thermal state of accreting white dwarfs experiencing 
classical nova explosions, Townsley and Bildsten (2004) have demonstrated that there 
exists a population of massive white dwarfs (M > I.2M0) accreting at rates as low as 
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«i 10~^^ MQyr~^ , with central temperatures Tc of ~ 5 x 10^ K. The combination of cold 
white dwarfs with low accretion rates leads to novae runaways with massive envelopes and 
exceptionally high peak temperatures. We show here that for a significant part of the 
relevant parameter space {Tc, Mwd and M) the peak temperature at runaway, Tpeofe is well 
above the minimum value for CNO "Breakout" . 

In section II we briefly describe the hydro solver and the reaction network. Section 
III presents the main results of our models. We conclude with a discussion of the possible 
observational consequence of our results and of their consistency with previously published 
models. 

2. The hydrodynamic solver 

The accretion process and the runaway are studied with a ID Lagrangian implicit code 
that integrates the equations of momentum and energy conservation assuming spherical 
symmetry. The energy transfer includes a radiative term and a convcctivc term. The 
radiative component is a diffusive flux. The radiative diffusion coefficient is determined 
according to the Iglasias & Rogers opacity tables Iglasias & Rogers 1996 for temperatures 
above Log{T) = 3.75 Kelvin and according to the Alexander fit for lower temperatures. 
Electron conductive opacities are computed according to the Itoh fit Itoh et. al. 1983. The 
convective energy flux is computed according to the mixing length theory Spiegel 1963. 
Within a convective zone the matter is mixed using a diffusion coefficient that takes into 
account both the local convective velocity and the mixing length (which is taken to be two 
pressure scale heights): 
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where Dc is the diffusion coefficient, Vc is the ID convective velocity, derived from the 
mixing length theory and I is the effective mixing length. 

In the present survey we ignore element settlement by diffusion and undershoot mixing. 
The equation of state employed to derive the energy and pressure for given temperature 
density and chemical abundances is an ideal gas for the ions and involves a table fit for 
the numerical integration of the Fermi integrals for the electrons (and positrons). The 
rate of nuclear energy production and the abundance changes are calculated with a 160 
nuclei network. The proton-proton reaction rate in the PeP channel take into account the 
extra electron density term due to the electron capture Bahcall & May 1969. The chemical 
elements included in the network were chosen in a proper way to provide an accurate 
history of the proton-rich composition of the accreted envelope. 

Accreted material is added each time step to the outermost zone as dictated by the 
accretion rate and the size of the time step. The outermost zone in the numerical grid is 
divided into two grid zones whenever the mass of the zone is greater than a specified value 
(taken to be twice the specified mass resolution in the envelope). The mass is added with 
the instantaneous thermodynamic properties of the outermost grid zone. The same ID 
solver was used in order to produce the initial models for all our 2D novae studies Glasner 
& Livne 1995, Glasner et al. 1997, Glasner et al. 2006. More details about the solver and 
its abilities will be given in a forthcoming paper. 

3. The main results 

In this study we follow the accretion process and the ensuing runaway for cold massive 
white dwarfs accreting at very low rates. Under such conditions the peak temperature 
achieved in the runaway is expected to be quite high. In order to trace the nucleosynthesis 
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history, we solve the hydrodynamic equations coupled with a nuclear reaction network that 
includes 160 nuclei up to ^^Ni. As we noted in the introduction, our interest in this poorly 
explored region of the nova parameter space was motivated by the results published by 
Townsley & Bildsten 2004. In our survey, we include white dwarf cores in the mass range 
1.2 Mq to 1.4 Mq accreting matter of solar composition Anders & Grevesse 1989 at rates 
ranging from 10~^^MQyr^^ to l{)^^^MQyr~^ . The initial models are cold white dwarf cores 
in hydrostatic and thermal equilibrium characterized by their central temperature (Tc) with 
values ranging from 7; = 4 x 10^ K to Tc = 2 x 10^ K. 

The combination of very low central temperatures Tc and very slow accretion rates, 
which demands the use of a greatly expanded nuclear reaction network, is the unique feature 
of this survey. In this short initial report we therefore focus on the main outcome: the 
extremely high peak temperatures achieved in the runaway and the induced nucleosynthesis. 
The main features of the time history of the total burning rate, the bolometric luminosity, 
and the peak temperature at the base of the envelope are quite similar for all studies 
models. 

Prior to the runaway, some of the considered models accrete exceptionally high 
envelope masses 6Mignite- The ensuing thermonuclear runaways are characterized by very 
high peak temperatures, Tpeak- In Figure 1 we present the relation between Tpeak and Tc 
for white dwarf masses and accretion rates (M) that are relevant to this work; in Figure 2 
we show the dependence of SMignUe on the same variables. The common feature for all the 
considered models is the increase of Tp^ak with decreasing Tci Tp^ak is greater than T^ru 
only for very cool white dwarfs (Tc < 5 — 7 x 10^). 

The high temperature and high density at the bottom of the accreted envelope and 
the comparatively long timescale for which these conditions prevail, defined as Stbreak, allow 
a significant epoch of explosive hydrogen burning. In Table I we present the main input 
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parameters of our models - Mwd, and M - together with the major hydrodynamic 
outcome quantities that are relevant to the burning: SMignUe , ^peofe a-nd Sthreak ■ 

As expected, hydrogen burning under such extreme conditions yields "breakout" 
from the conventional "hot carbon-nitrogen-oxygen" cycle. During most of the time that 
the temperature stays above Tcru the envelope is fully convective. Therefore, the nuclei 
produced in the burning are spread homogeneously over the entire envelope. Once the 
shell temperature falls below Tcru the overall abundance of heavy nuclei is not expected to 
change during the later stages of the runaway. Wc selected one representative model from 
Table I (sl35u for M=1.35 Mq) and present its initial and final abundance patterns as a 
function of mass number (A) in Figure 3. We note that both the mass of heavy elements 
ejected and their distribution depend on the details of the mass ejection mechanism; we 
will present a study of the ejection phase in a forthcoming paper. 



4. Discussion 

There is a well established consensus in the nova research community that the relevant 
parameters that define the thermonuclear runaway are the white dwarf mass (Mwd), the 
central temperature (Tc), and the accretion rate (M), together with the composition of 
both the cold white dwarf and the accreted matter. This parameter space has been well 
explored by many researchers. Most of the previous studies were concerned with the range 
of central temperature above approximately Tc = 10^ K. But, in a recent study, Townsley 
and Bildsten (2004) demonstrate that for very slow accretors there exists a steady state 
that determines a relation between M and Tc. For M less than 10~^°MQyr~^, they find 
that Tc according is much smaller than 10^ K. We therefore focused our study on massive 
white dwarfs with Tc much below Tc=10^ K and very slow accretion rates. 
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Since the expected temperatures at the runaway are above Tcru, the nuclear energetics 
and associated nucleosynthesis in our models were calculated with the use of an extended 
nuclear reaction network that included 160 nuclei up to ^'^Ni. Our results are generally 
consistent with those of the 1.3 Momodel of Nariai et al. (1980), as is evident from our 
Table I. For Tc=10^ K, both the accreted mass we obtain and the peak temperature Tpeak 
achieved during the runaway arc consistent with the 1.35 M0model of Jose and Hernanz 
(1998) and with the 1.25 Moand 1.40 Momodels of Yaron et al. (2005). For the conditions 
of central temperature and accretion rate we took from the results of Townsley and Bildsten 
(2004), the ignition mass we obtained is fully consistent with their stated results (compare 
our model sl25y in Table I with their Figure 8). 

Considering the dynamical aspects of our survey, we find a very steep dependence 
of the accreted mass and peak temperature at runaway on the central temperature Tc 
below Tc=10^ K. In contrast, the dependence in the more commonly studied region above 
Tc=10^ K is very shallow. This behavior can be understood on the basis of the well known 
dependence of the burning rates on the temperature during the accretion stage. For the 
accreted matter of assumed solar composition, the slope of the logarithmic derivative of the 
burning rate changes in the vicinity of a temperature 10^ from 4 (PP cycle) to about 18 
(CNO cycle). At central temperatures = 4 x 10^ K, relevant to cataclysmic variables 
below the period gap Townsley & Bildsten 2004, the accreted envelope mass required to 
initiate runaway can be 4 to 5 times greater than for Tc — 10^ K (Table I). This result 
confirms that once the core of the white dwarf is sufficiently cool, massive hydrogen-rich 
accreted envelopes can exist at runaway even for very massive white dwarfs. Therefore, 
observational evidence for nova outbursts with massive accreted envelopes in the range of a 
few times 10~^ MQcannot exclude very massive white dwarfs as the underlying core. 

Most significant are the results we obtain for the nucleosynthesis occurring in our 
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models. The high temperatures, much above T^rit, that can prevail for thousands of seconds 
lead to a very interesting path of nuclear burning up to the iron peak region. In order to 
identify the pattern of enrichment of iron group elements, we divided the elements included 
in our models into five groups: (1) Hydrogen; (2) Hehum; (3) all the CNO cycle elements 
up to fluorine; (4) the elements from neon to calcium (the intermediate mass elements); and 
(5) the iron group elements from titanium to nickel. For all the models reported in Table 
I, the peak temperature in the runaway exceeds Tcrit and produces large concentrations 
of iron group elements. The concentrations of freshly synthesized iron group elements 
(relative to the original solar amounts) range from an increase by a factor of 1.75 for the 
less massive models (model sl25y with M=1.25 Mq) to an increase by a factor of 10.0 for 
the more massive models (model sl35y with M=1.35 Mq); this can be seen from the results 
presented in table II. Both CNO elements (from breakout) and intermediate mass elements 
are seen to be transformed to iron group elements. The outcome is a net decrease in the 
mass of hydrogen and CNO elements. 

As one might expect, for the higher mass cases (higher peak temperature) there is even 
a small decrease in the fraction of intermediate mass elements and a significant production 
of iron group elements (see e.g. Figure 3 and Table II). For lower mass white dwarfs (lower 
peak temperature) there is a shght increase in the mass fraction of intermediate mass 
elements and a smaller level of production of iron group elements (see e.g. Table II). 

In order to make certain that we indeed achieved "breakout" from the conventional 
"hot carbon-nitrogen-oxygen" cycle we made a definite test. A model similar to sl35u 
(M=1.35 Mq), for which the accreted matter was assumed to include solar concentrations 
of nuclei only up to fluorine, was evolved. The energetics of this model are similar to those 
of model sl35u (Table I). With regard to the nucleosynthesis, the final concentration of iron 
group elements was determined to be comparable to that of model sl35u. There is only a 
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small concentration of intermediate mass elements in the final stages, which confirms that 
indeed CNO nuclei were burned all the way to iron. 

Our main conclusion concerning nucleosynthesis is that due to "breakout" from the 
conventional "hot carbon- nitrogen-oxygen" cycle, the heavy element abundances in the 
relevant novae ejecta can show a much wider variety than what is possible with the 
conventional mixing mechanisms. Some rare cases could be enriched by breakout alone. 
Future work will therefore concentrate on the observational predictions that can be made 
about the time history of the ejected mass. In order to test this hypothesis we will 
also consider the interpretation of observed abundance patterns in such nova systems as 
Nova Cygni 1992 and Nova Her 1991, in light of the possibihty that they may well have 
experienced nuclear burning at temperatures approximating or above breakout leading to 
increased nucleosynthesis of intermediate mass and heavy elements. In order to improve 
the accuracy of such predictions, we intend to follow the abundance history of the ejected 
matter all through the phase of nuclear burning and mass ejection. 
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Table 1. The main parameters for all models that pass T^ru 



[ MODEL] 


MwD [Mq] 


Tc [IQO K] 


M [Meyr-i] 


SMignite [Mq] 


Tpeafe [10« K] 


Stbreak [sec] 


sl35u 


1.35 


4.00 


10-11 


6.38 X 10-'"^ 


5.14 


14500. 


ml35k 


1.35 


5.88 


10-11 


3.90 X 10-'"^ 


4.57 


12000. 


ml35j 


1.35 


7.65 


10-11 


2.70 X lO-'"^ 


4.23 


9500. 


ml35i 


1.35 


9.03 


10-11 


2.06 X 10-'^ 


4.00 


7000. 


ml35r 


1.35 


4.00 


10-10 


3.40 X 10-5 


4.42 


11000. 


ml35q 


1.35 


5.88 


10-10 


2.63 X 10-5 


4.16 


9000. 


nil35p 


1.35 


7.65 


10-10 


2.01 X 10-5 


3.92 


7000. 


nil35o 


1.35 


9.03 


10-10 


1.61 X 10-5 


3.51 


1000. 


ml3ek 


1.30 


4.00 


10-11 


9.62 X 10-5 


4.29 


7000. 


mlSed 


1.30 


6.00 


10-11 


5.41 X 10-5 


3.82 


4250. 


sl25y 


1.25 


4.00 


10-11 


13.02 X 10-5 


3.68 


3000. 
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Table 2. The changes in mass abundances grouped according to mass number for the 

models sl35u and sl25y 



jobname 


group of elements 


absolute change 


relative change 


sl35u 


AX[H] 


-0.0740 


-0.10 




AX [He] 


+0.0660 


+0.24 




AX[C-F] 


-0.0023 


-0.23 




AX[Ne - Ca] 


-0.0032 


-0.90 




AX[Ti-Ni] 


+0.0135 


+10.0 


sl25y 


AX[H] 


-0.0377 


-0.05 




AX [He] 


+0.0365 


+0.13 




AX[C - F] 


-0.0007 


-0.07 




AX[Ne - Ca] 


+0.0009 


+0.25 




AX[Ti-Ni] 


+0.0010 


+0.75 
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Fig. 1. — The peak temperature achieved in the runaway as a function of the central 
temperature Tc of the white dwarf. 
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Fig. 2. — The accreted mass at ignition as a function of the central temperature of the 
white dwarf. 




Fig. 3. — The abundances of the elements in the reaction network as a function of the mass 
number A. Presented are both the initial (solar) abundances and the 'final' abundances when 
the peak temperature has fallen below Tcru (Model sl35u, M=1.35 Mq). 



